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Abstract 
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' We examine the phenomenological implications of a ~ 200 GeV leptoquark in light 

m ■ 

£T) • of the recent excess of events at HERA. Given the relative predictions of events rates in 

m : 

e + p versus e~p, we demonstrate that classes of leptoquarks may be excluded, includ- 
ing those contained in Eq GUT models. It is shown that future studies with polarized 

Oh, beams at HERA could reveal the chirality of the leptoquark fermionic coupling and 

Oh! , 

that given sufficient luminosity in each e L R channel the leptoquark quantum numbers 

43 ■ 

could be determined. The implications of 200 — 220 GeV leptoquarks at the Tevatron 
are examined. While present Tevatron data most likely excludes vector leptoquarks 
and leptogluons in this mass region, it does allow for scalar leptoquarks. We find that 
while leptoquarks have little influence on Drell-Yan production, further studies at the 
Main Injector may be possible in the single production channel provided the Yukawa 
couplings are sufficiently large. We investigate precision electroweak measurements as 
well as the process e + e~ — ► qq at LEP II and find they provide no further restrictions 
on these leptoquark models. We then ascertain that cross section and polarization 
asymmetry measurements at the NLC provide the only direct mechanism to determine 
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the leptoquark's electroweak quantum numbers. The single production of leptoquarks 
in 7e collisions by both the backscattered laser and Weisacker- Williams techniques 
at the NLC is also discussed. Finally, we demonstrate that we can obtain successful 
coupling constant unification in models with leptoquarks, both with or without super- 
symmetry. The super symmetric case requires the GUT group to be larger than SU(5) 
such as flipped SU(5) x U(l)x- 
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1 Introduction 



The apparent symmetry between the quark and lepton generations is a mysterious occurrence 
within the Standard Model (SM) and has inspired many theories which go beyond the SM 
to relate them at a more fundamental level. As a result many of these models naturally 
contain leptoquarks, or particles that couple to a lepton-quark pair. Theories which fall in 
this category include, composite models with quark and lepton substructure^], the strong 
coupling version of the SM0, horizontal symmetry theories^, extended technicolor and 
grand unified theories (GUTs) based on the gauge groups SU (5)0, SO (10) with Pati-Salam 
S77(4) color symmetry ||, 577(15)117), an d superstring-inspired E% models||, [J. In all cases, 
the leptoquarks carry both baryon and lepton number and are color triplets under SU(3)o 
In models where baryon and lepton number are separately conserved, which includes most 
of the above cases, leptoquarks can be light (of order the electroweak scale) and still avoid 
conflicts with rapid proton decay. Their remaining properties, such as spin, weak iso-spin, 
electric charge, chirality of their fermionic couplings, and fermion number, depend on the 
structure of each specific model. If leptoquarks were to exist we would clearly need to 
determine these properties in order to ascertain their origin. 



An excess of events at large values of Q 2 have recently been reported Jl0| by both 



the HI and ZEUS collaborations at HERA in their neutral current Deep Inelastic Scattering 
(DIS) data. ZEUS has collected 20.1 pb _1 of integrated luminosity in e + p collisions and 

observes 5 events with Q 2 > 15, 000 GeV 2 with x > 0.45 and y > 0.25, where x and y are 
the usual DIS scattering variables, while expecting 2 events from the SM in this region. 
HI reports 7 events in the kinematic region m = \fxs > 180 GeV and y > 0.4, compared 

to a SM prediction of 1.83 ± 0.33 with 14.19 ± 0.32 pb -1 of integrated luminosity. Clearly, 
the statistical sample is too small at present to draw any conclusions and it is likely that 
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this excess is merely the result of a statistical fluctuation. Another possibility is that this 
discrepancy is the result of deviations from current parton distribution parameterizations at 
large x. This case, however, has been examined by the HI and ZEUS collaborations [|TU[ and is 
found to be unlikely. Also, such large modifications in the parton densities would most likely 
result in disagreement with the dijet data samples at the Tevatron[[TT|, [12]. It is also possible 



that this HERA data might signal the first hint of physics beyond the SM. Such an excess in 
event rate at large Q 2 is a classic signature for compositeness if the events show no specific 
kinematic structure. This scenario has recently been analyzed |13| in light of the HERA 



data, with the result that an eeqq contact interaction with a right-left helicity structure (in 
order to avoid the constraints arising from Atomic Parity Violation data discussed below) 
and a scale of ~ 3 TeV is consistent with the data. However, if instead, the events cluster 
in x while being isotropic in y, they would signal the production of a new particle. While 
the reconstructed values of the mass (m = y/xs) of such a hypothetical particle show some 
spread between the two experiments, they are consistent within the evaluated errors, yielding 
a central value in the approximate range 200 — 210 GeV. 

If this excess of events turns out to be the resonant production of a new particle, we 
will need to examine the possible classes of new lepton-hadron interactions which could give 
rise to such a signature. At present, there are three leading scenarios of this type: (i) models 
with leptogluons, (ii) explicit R-parity violating interactions in Supersymmetric theories, 
and (iii) leptoquarks. We now briefly discuss the first two cases, with the remainder of the 
paper being devoted to the third. 

Color octet partners of ordinary leptons are expected to exist in composite models 
with colored preons||14|. These particles, denoted as leptogluons, are fermions, carry lepton 



number and couple directly to a lepton-gluon pair with an undetermined strength. This 
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effective interaction may be written as 

Ce,f = | [\ L La^£l L + X R £l R a^e R ] G%£ + h.c. , (1) 

where G a is the gluon field strength tensor, A is the compositeness scale, L represents the 
lepton doublet under SU(2)l, and Xl,r parameterizes the unknown coupling. We note that 
this effective Lagrangian is non-renormalizable. If the chiral symmetry in these models is 
broken by QCD effects, the leptogluons are expected ||14|| to have masses or order a s A, and 



hence could be as light as a few hundred GeV for compositeness scales in the TeV range. 
Leptogluons in this mass range would clearly reveal themselves in high-Q 2 DIS at HERA. A 
~ 200 GeV leptogluon would also be copiously produced at hadron colliders. In fact, as we 
will see below, present Tevatron data most likely excludes leptogluons in this mass region. 

The second scenario for resonant new particle production listed above is that of 
super symmetric theories with explicit R-parity violating interactions. The most general 
gauge and supersymmetry invariant superpotential (with minimal field content) contains 
the terms 

XijkLiLjEk + X'ijfrLiQjDk + X'-j k UiDjDk , (2) 

where the first two terms violate lepton number (L), the third violates baryon (B) number, 
the A's are a priori unknown Yukawa coupling constants, and the i,j, k are generational 
indices where SU(2)l invariance demands that i ^ j. In the Minimal Supersymmetric Stan- 
dard Model (MSSM) a discrete symmetry matter parity (or R-parity) is applied to prohibit 
all these dimension four B and L violating operators. However, it is sufficient to ensure only 
that the B and L violating terms do not exist simultaneously in order to preserve nucleon 
stability. In fact, a discrete anomaly free Z% symmetry, denoted as baryon parity, naturally 
allows for the L violating operators, while forbidding the AB ^ operators (TTSJ. The phe- 
nomenology of these models is strikingly different than in the MSSM, as elementary couplings 



involving an odd number of supersymmetric particles now exist. This results in the possible 
single production of super-partners and an unstable Lightest Supersymmetric Particle. At 
HERA, the second L violating term in Eq. Q can mediate single squark production. This 
possibility was first considered in Ref. [16|, and later examined in detail in Ref. [17]. The 



relevant terms in the interaction Lagrangian for this case are 

C = -X' ljk [u{d k R e l L + (4)*(e!)V L ] + h.c. (3) 

The requirement of SU(2)l invariance combined with the fact that the positrons must be 
scattered off of valence quarks in order to account for the event excess at HERA, leaves 
us with only one possible scenario, the production of charm or stop squarks via the first 



term above. This case has been recently examined [18] in light of the data and will not be 
considered further here. However, we note that such singly produced squarks must also 
decay via their R-parity conserving interactions (e.g., q — > q + x°) a t competitive rates, and 
hence events with these signatures must also be observed. We will comment on some of the 
difficulties associated with GUT theories with R-parity violation below. 

We now examine the third candidate scenario above, the existence of 200 — 210 GeV 
lepto quarks, in detail. 

2 What is a Leptoquark? 

The interactions of leptoquarks can be described by an effective low-energy Lagrangian. 
The most general renormalizable SU(3)c x SU(2)l x U(l)y invariant leptoquark-fermion 



interactions can be classified by their fermion number, F = 3B + L, and take the form[]19 

C = C F= ^ 2 + C F=0 , (4) 



5 



with 

£f=-2 = {giLql^L + gmu c R e R )Si + giRd C R e RSi + g3Lq c L ir 2 r£ L S 3 

Hg2Ld c R l/L + g2Rq c L Ye R )V 2fl + g 2 LU c R YhV 2lI + h.c. , (5) 

£f=o = (h 2L u R £ L + h 2R q L ir 2 e R )R 2 + h 2L d R £ L R 2 + (h 1L q L Yh + hi R d R ^c R )Ui^ 

+h 1R u R n lx e R U llx + h ZL q L f^i L U^ + h.c. 

assuming the fermionic content of the SM. If the SM fermion content is augmented it is 
possible that new LQ interactions may arise. Here, qi and £l denote the SU(2)l quark 
and lepton doublets, respectively, while u R , d R and e R are the corresponding singlets. The 
indices of the leptoquark fields indicate the dimension of their SU{2)l representation. The 
subscripts of the coupling constants label the lepton's chirality. For simplicity, the color 
and generational indices have been suppressed. Since, in general, these couplings can be 
intergenerational, there is the possibility of large, tree-level flavor changing neutral currents 
and flavor universality violations. As discussed in the next section, this can be avoided by 
employing the constraint that a leptoquark couple only to a single generation. We see that 
the leptoquark fermionic couplings are baryon and lepton number conserving, hence avoiding 
the conventional problems associated with rapid proton decay. Note that the leptoquarks 
with fermion number (F) of —2 (S and V) couple to £q, while the F = leptoquarks (R 
and U) have £q couplings. Once this effective Lagrangian is specified, the gauge couplings 
of the leptoquarks are completely determined. Thus, only the strength of the leptoquark's 
fermionic Yukawa couplings remain unknown. For calculational purposes, these couplings 
are generally scaled to the electromagnetic coupling, 

Al„r = e\ LtR , (6) 

where \l,r generically represent the giL, R and hiL,R. We note that in extended technicolor 
theories, the leptoquarks (denoted there as P3) have couplings which are proportional to the 
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masses of the quark-lepton pair. In this case the qeP% coupling is clearly too small to be of 
interest to us here. 

The quantum numbers and structure of the quark-lepton couplings are summarized in 
Table 1 for each leptoquark species. For the couplings we list the helicity, relative strength, 
and fermion pairs with which a particular leptoquark may couple, using the convention 
LQ F= _ 2 — > £q and LQ F=0 — > £q. The leptoquark's branching fraction into charged leptons, 
Bi, is also given. The weak isospin structure is denoted by the brackets. Due to gauge 
invariance we would expect all the leptoquarks within a given SU(2)l representation to 
be degenerate apart from loop corrections. For future reference we have also listed the 
SU(5) representations of lowest dimension within which the leptoquark can be embedded. 
Note that we have assumed the conventional assignments of the SM fermions to be in the 
5 and 10 representations. Generally, only a subset of these possible leptoquark states are 
contained within a particular model. For example, the scalar Sh,r is the leptoquark present 
in superstring-inspired Eq theories. One exception is the GUT based on 577(15), which 
contains all 14 possible leptoquark states! 

Using Eq. 5, the total leptoquark tree-level decay widths are easily calculated to be 

r= E; Af = 0.386 ^q^^E^ GeV, Scalars 

Ej A? = 0.258 20() ^ e y J2i Af GeV , Vectors 

where we have scaled the widths to a 200 GeV leptoquark in our numerical evaluation, and 
the sum extends over all possible decay modes. These states are clearly very narrow and 

hence are long-lived, especially for the values of A that are consistent with the low-energy 
constraints discussed in the next section. As pointed out by Kunszt and Stirling, QCD 
corrections to these widths are very small . 



7 





Lepto quark 


SU(5) Rep 


Q 


Coupling 


B e 


Scalars 




F = -2 


•SlL 


5 


1/3 


Al (e+u), X L (vd) 


1/2 




Sir 


5 


1/3 


\ R (e+u) 


1 




Sir 


45 


4/3 
4/3 


/v il. V c u / 

-^/2A L (e+d) 


1 
1 






45 


1/3 
-2/3 


-Al (e+u), -Al (Pd) 
V2A L (Pu) 


1/2 



F = 


R2L 


45 


5/3 
2/3 


A L {e+u) 
Al (^«) 


1 





R2R 


45 


5/3 
2/3 


Ab (e+u) 
-Xr (e+d) 


1 
1 




R2L 


10/15 


2/3 


Al (e+d) 
Al {yd) 


1 




Vectors 




F = -2 


v 2L 


24 


4/3 
k 1/3 


Al (e+d) 
A L {yd) 


1 






V 2 R 


24 


4/3 
k V3 


Xr (e+d) 
Ar (e+u) 


1 
1 




v 2L 


10/15 


1/3 
-2/3 


Al (e+u) 
Al (^m) 


1 




F = 


u 1L 


10 


2/3 


Al (e+d), Al (Pu) 


1/2 




UlR 


10 


2/3 


A_r (e+d") 


1 




UlR 


75 


5/3 
5/3 


Xr (e+u) 
V2A L (e+u) 


1 
1 




U 3 L 


40 


2/3 


-Al (e+d), Al (^u) 


1/2 






8 


-1/3 


V2Al (i/d) 






Table 1: Quantum numbers and fermionic coupling of the leptoquark states. No distinction 
is made between the representation and its conjugate. 



3 Low Energy Constraints 



As mentioned above, low-energy data places strong restrictions on the leptoquark Yukawa 
couplings. In this section we summarize the most relevant of these constraints. 

As is well known, for a leptoquark to be sufficiently light for it to be of phenomeno- 
logical interest at existing or planned colliders, it must have essentially chiral couplings to 
fermions. Here, we give two examples which demonstrate this conclusion, (i) Consider a 
first generation leptoquark coupling to ue^ with strength \ R , and to dv L with strength \ L . 
A Fiertz transformation then yields the interaction 

^ = u Rd L e R is L , (8) 

2m LQ 

which gives a large contribution to the decay 7r + — > e + u e . Comparison with current data 
results [pi], |22| in the bound m L Q > 200\\ L \ R \ 1 / 2 TeV, implying that at least one if not 



both of these Yukawa couplings must be small, if leptoquarks are to be light. If one of 
these couplings is sufficiently large in order to induce leptoquark-fermion interactions at an 
interesting level, then the coupling with the other handedness must essentially vanish, i.e., 
the couplings are chiral. (ii) One can also examine the leptoquark's contribution to the 
g — 2 of the muon arising from a one loop penguin diagram involving a light quark and 
a leptoquark. In this case if couplings of both helicities are present, current data places 
the constraint^]]] Mlq > IOOOAlAj? TeV, and again we see that these couplings must be 
essentially chiral. 

Even if the leptoquark-fermion couplings are chiral, strong constraints on their mag- 
nitude still arise from their potential contributions to a wide class of flavor changing neutral 
currents and related phenomena. An exhaustive study of this class of transitions has been 
performed by Davidson, Bailey and Campbell||21||, and hence we will not repeat this type of 
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investigation here. The results of this study show that most of the stronger bounds can be 
satisfied if the requirement that a given leptoquark couple to only one generation is imposed. 
This results in the nomenclature of first, second and third generation leptoquarks found in 
the literature. Note that this restriction is more severe than the simple requirement of family 
number conservation at a single vertex. This is easily illustrated by examining the process 
K — > fie, which would receive a large tree-level contribution if leptoquarks were allowed to 
simultaneously couple to both the first and second families. 

After the constraints of chiral and single generation couplings are imposed, there are 
two important remaining low-energy constraints arising from Atomic Parity Violation (APV) 
and the universality testing decay tt ev. Comparable but somewhat weaker bounds also 
follow from quark-lepton universality. These additional restrictions have been examined by 
Leurer [p2| and also by Davidson et a/- ||21|| for both cases of scalar and vector leptoquarks. We 



summarize these results in Table 2 (assuming that the LQ is not responsible for the small 
difference between the SM expectations for the APV 'weak charge' and what is obtained 
experimentally). These values have now been updated to include the recent results of Wood 
et al.on APV in Cesium[]23| which are in good agreement with SM predictions |24| yielding 



^Qw — 1-09 ± 0.93. Note that these bounds are far from trivial. For example, we see from 

this Table that a spin-0, i?2L leptoquark with a mass of 200 GeV must have A < 0.22! As 
we will see below this is not far from the value suggested by the excess of events at HERA. 

Other types of experiments give slightly weaker bounds on the Yukawa couplings for 
fixed leptoquark mass. For example, it is well-known that precision measurements in deep 
inelastic neutrino scattering are sensitive to new particle exchanges. Using the latest CCFR 

results^] we obtain the constraint \l,r ^ 0.4 — 0.6 for 200 GeV leptoquarks. NuTeVfl2"o| 
may be able to improve this reach by a factor of two. Older results, such as that from the 



SLAC polarized electron-Deuteron scattering experiment f26|1 , are only sensitive to \{l,r of 
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Table 2: Combined limits on the ratio m/X in GeV, where m is the leptoquark mass, for the 
leptoquarks multiplets from data on Atomic Parity Violation and the decay ir — > ev. 



order unity or greater. 



4 Leptoquarks at HERA 

Clearly, ep collisions are especially well suited for leptoquark production and offer striking 



signals|19|, |27|, |28|. Direct production contributes to DIS in either the neutral or charged 
current channel, through an s-channel resonance e ± ^ — > LQ with the subsequent decay to 
either e ± (q> or v e (q> with a fixed branching fraction, depending on the leptoquark species. 
The final state with v e manifests itself as missing energy. Clearly, this s-channel exchange 
would yield distinctive, and due to the size of the width, narrow peaks in the x-distributions 
at x = m 2 /s. These peaks, however, are smeared by the detector resolution as well as QCD 



and QED radiative effects ]29|, pOj . Additional smaller contributions are generated from u- 
channel leptoquark exchange. The fermion number of the leptoquark dictates whether it 
will contribute via s- or u-channel exchange in e~ versus e + scattering off of valence (q) or 
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sea (q) partons. For example, the F = — 2 leptoquarks (S and V) mediate DIS through the 
s-channel (u-channel) in e + q (e + q) collisions, while the F = states are exchanged in the 
u-channel (s-channel) in e + q (e + q) collisions. In principle, this can be used to separate the 
production of F — —2 from F = leptoquarks from cross section measurements alone. 

We first examine the expected event yield in the e ± j channel at HERA for the pro- 
duction of the various leptoquark species. We concentrate on the case of scalar leptoquarks, 
as ~ 200 GeV vector leptoquarks are most likely excluded by Tevatron data as shown in the 
next section. The cross section is dominated by the s-channel resonance, and the narrow 
width of these states as seen as Eq. [7| justifies the use of the narrow width approximation. 
The differential cross section for leptoquark production can then be written as 

dajep ~^LQ^ ej) n 2 a~ 2 H 2 2 J 1, Scalar! 
= A ( q>(m /s,Q )B e < }. (9) 

d V s { Q{l-y) 2 , Vector 

For the case of scalar leptoquarks, we see that the production is isotropic in y, whereas the 
electroweak DIS background has a 1/y 2 behavior. In obtaining the total cross section for 
scalar leptoquarks, we use the y-averaged parton densities 

(»(«,«»)> = Jg(x ;^ . do) 



with the range 0.25 < y < 1.0, and employ the MRSA' distributions |3T[]. The resulting 



expected excess of events from scalar leptoquark production are displayed in Table for 
e + p and e~p collisions, scaled to 20 pb -1 and 1 pb _1 of integrated luminosity, respectively. 

Here we assume m = 200 GeV and for purposes of demonstration we take A = 0.1. The 
numbers in brackets represent the corresponding event rate in the faj channel. Note that 
only 2 leptoquark species can contribute to charged current DIS assuming only the SM 
fermion content as shown in Eq. 5. We stress that the relative magnitudes of these event 
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Table 3: Number of events for scalar leptoquark production in e + p and e p collisions, scaled 
to 20 pb -1 and 1 pb -1 of integrated luminosity, respectively, assuming m = 200 GeV, A = 
O.le, and taking 0.25 < y < 1.0. The numbers in brackets indicate the corresponding 
expected event yield in the charged current channel. 

rates are fixed by the contributing parton densities and the luminosity. As expected, the 
Si leptoquarks have significantly larger cross sections in electron (rather than positron) 
collisions, since the valence quark distributions contribute in this case. Thus, in order to 
account for the HERA data, the F = —2 scalar leptoquarks would also yield an excess of 
events in the ~ 1 pb -1 of e~p data. For example, if the leptoquark contained in E 6 theories, 
Sil, were to account for the observed e + j excess, then it would also yield roughly 600 excess 
events in 1 pb -1 of e~p data! Thus, unless there are significant event excesses hiding in the 
HI and ZEUS e~p data (which is not yet completely analyzed), we may exclude F = — 2 
S-type leptoquarks as the source of the HERA events. In contrast, F = scalar leptoquark 
production is suppressed in electron collisions. We see that in the case of R 2 l, the predicted 
number of events in the e + j channel with the assumed coupling strength of A ~ 0.1 is 
consistent with the data, while for R 2 l,r the coupling would have to be somewhat smaller 

with A ~ 0.03 — 0.04 neglecting the potentially large QCD corrections. 

If the leptoquark signature is verified by future data taking at HERA, it will be 
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mandatory to determine its couplings. Clearly, the best method of accomplishing this at 
HERA is to use both e ± p collisions and to take advantage of possible beam polarization [19, 



|28fl . (It is expected that polarization levels of P ~ 50% may be achievable at HERA in the 
future.) Table § displays the total number of expected events, scaled to 100 pb~ x , assuming 
100% beam polarization for ef R p collisions for a 200 GeV scalar leptoquark of each type, 
subject to the cuts 0.4 < y < 1.0 and M e j = 200 ± 20 GeV to remove the SM background. 
In these calculations the full deep inelastic scattering amplitudes, including the exchanges 
in all channels, have been used. The numerical results justify our earlier use of the narrow 
width approximation. The results have also been smeared by a 5% mass resolution. It is 
clear from the Table that knowledge of the ratio of cross sections which are essentially free 
of QCD corrections in the four channels will allow the leptoquark quantum numbers to be 



determined if sufficient statistics are available. For fixed values of A, it has been shown [20|. |32|| 
that the QCD corrections to the production of leptoquarks off of the valence partons can be 
as large as +25%, but are somewhat smaller for the case of production off of sea quarks. 

5 Signatures at the Tevatron 

Leptoquarks may reveal themselves in several reactions at hadron colliders. They may be 
observed directly via pair or single production mechanisms or they may indirectly influence 
the lepton pair invariant mass spectrum in Drell-Yan processes. We examine each of these 
in this section. 

Since leptoquarks are color triplet particles, their pair production |33|, |34j] proceeds 



through gluon fusion or quark annihilation and is essentially independent of the Yukawa 
coupling, A. There is a potential contribution of order A 2 via the reaction qq — > LQLQ with 
t-channel lepton exchange, however this contribution is negligible for the size of Yukawa 



14 



Leptoquark 


NE 




Nt 




SM background 


51.7 


28.7 


9.98 


20.0 


SlL 


121. 


28.7 


9.98 


20.4 


Sir 


51.7 


167. 


10.8 


20.0 


Sir 


51.7 


63.0 


11.5 


20.0 


Ssl 


190. 


28.7 


9.98 


23.5 


R2L 


52.4 


28.7 


9.98 


158. 


R2R 


51.7 


29.4 


148. 


20.0 


R2L 


53.2 


28.7 


9.98 


54.4 



Table 4: Number of events per 100 1 pb for each electron charge and state of polarization for 
a 200 GeV scalar leptoquark at HERA assuming 0.4 < y < 1, A = 0.1, and M ej = 200 ± 20 
GeV. These results have been smeared with a detector resolution of 5% in M e j. 

couplings of relevance here, A ~ C?(10 _1 e). The pair production of scalar leptoquarks thus 
mimics that of squarks. The number of events expected at the Tevatron, scaled to 100 
pb -1 of integrated luminosity, for the pair production of one generation of a single type of 
scalar leptoquark is displayed as the solid curve in Fig. [1]. Here, we have employed the 
MRSA' parton distributions []3l], and omitted the i^-factor which has been calculated |35| for 
leptoquark production to be K gg = 1 + 2a s 7r/3 and K qq = 1 — a s 7r/6 for the gluon fusion 
and quark annihilation subprocesses, respectively. For a 200 GeV scalar leptoquark, this 
yields an enhancement in the cross section by a factor of 1.16 giving a = 0.117 pb assuming 
/i 2 = s. If instead fi 2 = m 2 is chosen, a larger cross section will result since we always have 
s > Am 2 . We note that the cross section falls rapidly, dropping by a factor of 1.32 between 
m = 200 and 210 GeV with a(m = 210 GeV) = 0.089 pb (including the K-factor). We note 
that a complete NLO calculation of scalar LQ pair production has now been completed by 
Kramer et al. which essentially reproduces the results obtained using the K-factor approach 
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at a scale of fi 2 = m? p6[ . 



The signatures for leptoquark pair production are 2 jets accompanied by either i + l , 
fe>T or i>T i with a pair of jet+£ ±,t} invariant masses being equal to the mass of the 
leptoquark. DO has searched for the dijet with two or single charged lepton topologies in the 
electron and muon channels, and CDF has searched in the dilepton+dijet case for all three 

are m,( e ) > 



lepton generations. For each generational coupling the most stringent bounds |p7 
175(147) GeV from DO, m (p) > 180(140) GeV and m (r) > 99 GeV from CDF, with B e = 
1(0.5). DO sets a 95% C.L. limit on the pair production cross section of scalar leptoquarks 
decaying into e + e~jj of a <^ 0.25 pb _1 . This DO bound is based on the observation of 3 
events with a Monte Carlo background estimate of 2.85 ± 1.08 events arising from Drell-Yan, 
tt — > ££, Z — > tt — > tt, WW -»■ it and fakes from QCD. Clearly the leptoquark SU(2) L 
representations which contribute to the cross section significantly more than that of a single 
leptoquark, e.g., R^r, will have a more difficult time satisfying these bounds. We stress 
again that these constraints are independent of the value of A. 

CDF has yet to present results from a search for first generation leptoquarks from 
runs 1A and IB, which contain a combined data sample of approximately llOpfe -1 . Based 
on our cross section above, we would expect CDF to observe at most 1 signal event (taking 
Be = 1) depending exactly on the leptoquark mass (e.g., 200 versus 210 GeV) and the 
details of their selection criterion, plus an unknown amount of background. Until all of these 
numbers become available we cannot attempt to combine the CDF and DO results to obtain 
a stronger mass bound without speculating upon the details of the CDF data and a thorough 
understanding of the common systematics of the two experiments. However, we would not 
expect a combined DO/CDF bound to significantly exceed 200 GeV. 

In order to compute the pair production cross section for vector leptoquarks (V) we 
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need to determine both the trilinear gVV and quartic ggVV couplings. In any realistic 
model that contains fundamental vector leptoquarks, they will be the gauge bosons of some 
extended gauge group. Hence gauge invariance will completely specify the gVV and ggVV 
couplings in such a manner as to guarantee that the subprocess cross section obeys tree-level 
unitarity, as is the hallmark of all gauge theories. However, vector leptoquarks could be a 
low-energy manifestation of a more fundamental theory at a higher scale, and they could 
be composite. In this case various anomalous gVV and ggVV couplings could be present, 
one of which can be described by a chromo-magnetic moment, k. This term represents the 
only dimension 4 anomalous coupling which conserves CP. In gauge theories k takes on the 
value of unity. The cross sections for vector leptoquarks with k = 1 have been computed in 
Refs. |38], In this case the total event rate, scaled again to 100 pb -1 of integrated 
luminosity, for VV production at the Tevatron is given by the dashed curve in Fig. [I]. Vector 
leptoquarks have the same signatures as discussed above for the scalar case, but with slight 
detailed variations in the production angular distribution due to the fact they are spin-1 
particles. A reasonable estimate of the search reach can be obtained by employing the DO 
bound on the cross section for e + e~jj events from scalar leptoquark production. We estimate 
that this procedure yields the constraint on first generation vector leptoquarks with k = 1 
of m ^ 290 GeV, placing this case out of the kinematic reach of HERA |f40|| . Clearly, the 
experiments themselves need to perform a detailed analysis in order to confirm this estimate. 



The results in the more general case|^, 38 1 of k ^ are displayed in Fig. [| Here 



the separate qq, gg, and the total cross sections for vector leptoquark pair production with 
m = 200 GeV are shown as a function of k. We see that the cross section varies significantly 
with k yielding larger or smaller values than the results given above for k = 1. In the worst 
case, the total cross section reaches its lowest value of <Ji ow m 0.6 pb around k ~ —0.45. 
This value of k is, of course, much larger than one would expect in any realistic model but 
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is considered here for generality. However, this value of o\ ow is significantly larger than the 



DO bound|37| of 0.25 pb and would again exclude 200 GeV vector leptoquarks for this n 
with Bi = 1. For all other values of k, the cross section is comfortably large enough to 
be prohibited by DO. We remind the reader that we have neglected the K-factors in this 
analysis; their inclusion will only strengthen our conclusions. 

Here, we also consider the case of pair production of leptogluons at the Tevatron. This 
process was considered in Refs. |4l| , f42f and is also mediated by gg fusion and qq annihilation, 
similar to the production of any heavy colored fermion. However, since leptogluons are color 
octets, there is an enhanced color structure in this case compared to, e.g., top-quark pair 
production. The event rate, scaled to 100 pb -1 , is given by the dotted curve in Fig. [I]. We 
see that this cross section is larger than that for both scalar and vector leptoquarks, and is 
roughly six times the top pair cross section. The DO cross section bounds on e + e~jj events 
would clearly exclude 200 GeV leptogluons and could naively place the constraint mic <L 325 
GeV. 

Leptoquarks can also be produced singly at hadron colliders. The parton level sub- 
process responsible for single production is qg — > LQ + £, where I is a charged lepton or 
neutrino depending on the type of leptoquark. The diagrams for this process contain the 
QCD strong coupling at one vertex and the leptoquark Yukawa coupling at the other. The 
subprocess differential cross section for scalar leptoquarks is given by|53|, Q 



da iia.aX 2 



dt 3s 5 



s + i-m 2 + i(i + m 2 ) + i(2m 2 - s) 



(t — m 2 ) 2 s(t — m 2 ) 



(11) 



where i = |(m 2 — s)(l — cos#) with 9 being the quark-lepton scattering angle. For com- 
pleteness we give the corresponding single production cross section for vector leptoquarks in 
the Appendix. Compared to pair production, this mechanism has the advantage of a larger 
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Figure 1: Pair production cross sections at the Tevatron for scalar and vector leptoquarks, 
as well as leptogluons, corresponding to the solid, dashed, and dotted curves, respectively. 
Here, the gg and qq contributions are summed. 
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Figure 2: k dependence of the qq, gg, and total vector leptoquark pair production cross 
sections at the Tevatron, represented by the dotted, dashed, and solid curves, respectively. 
The leptoquark mass is taken to be 200 GeV. 



20 



amount of available phase space, but has the disadvantage in that it is directly propor- 
tional to the small Yukawa coupling. The total cross section at the Tevatron Main Injector 
and/or TeV33(now taking y/s = 2 TeV) in the case of scalar leptoquarks with A/e = 1 for 
both gu + gu and gd + gd fusion is presented in Fig. [| Note that for a pp collider, the 
gq and gq cross sections are equal. Here, we have again omitted the K-factor (given by 
K = 1 + 3« s 7r/4 for single production|[]|), which enhances the cross section by a factor 
of ~ 1.25. For a 200 GeV scalar leptoquark with coupling strength of A/e = 0.1 and in- 
cluding the K-factor, our results show that we obtain approximately ~ 39 , 87 events from 
gd + gd , gu + gu fusion, respectively, would be obtained with 10 fb" 1 of integrated luminosity 
at the Main Injector/TeV33. This event rate should be sufficient to provide a very rough 
determination of the value of Yukawa coupling A. The signatures for this production mode 
are jet +£ + £~ , +£ ± $t ,+ i>T and at least in the first case, should be easily detectable. 

Indirect signals for leptoquarks may be observed in Drell-Yan production ]33J]. In 
addition to the usual s-channel 7 and Z exchange in the SM, leptoquarks may also contribute 
to qq — ► e + e~ in t- or u-channel exchange, with the specific channel depending on the 
leptoquark type, via the Yukawa coupling. This new exchange will modify not only the 
invariant mass distribution, but also the angular distribution of the lepton pair. The parton 
level differential cross section for this process can be written as 



where M represents the invariant mass of the lepton pair, z = cos 9* with 9* being the 
center-of-mass scattering angle, K is the usual QCD correction factor with a s evaluated at 
the scale M, and the sum extends over the appropriate partons. In the SM all quark flavors, 
in principle, contribute to this process, whereas, in the case of leptoquark exchange, only 



da 



Knot 2 



Y^\- F q G q + F q G q]y 



(12) 



dMdydz 



6M 3 
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Figure 3: Single production cross sections at the Tevatron for scalar leptoquarks as a func- 
tion of mass from gu{u) and gd(d) fusion, corresponding to the dotted and dashed curves, 
respectively. Here, A/e = 1, and the K-factor has been omitted. 
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one or two quark flavors contribute. The parton density factors are given by 



= x a x b [q{x a )q{x b ) ± q{x b )q{x a )} 



(13) 



with x a>b = (M/y / s)e ±J/ as usual, and represent the even and odd kinematic functions, 
which are given in the Appendix for both the SM and leptoquark contributions. 

The invariant mass distribution is obtained by integrating the differential cross section 
over the regions — Y < y <Y and — Z < z < Z where 



Y = min[y max ,-ln{M/yG)] 
Z = mm[tanh(y — \y\), 1] 



(14) 



where y ma x represents the rapidity coverage of the detector or of the applied cuts. To calculate 
the forward-backward asymmetry, Eq. ^ is first integrated over the forward (z > 0) and 
backward {z < 0) regions separately (subject to \z\ < \Z\) and then over y; the difference of 
the forward and backward cross sections divided by their sum then gives Afb- Explicitly, 
we define 



da + 
dM 



da 
~dM 



f dy+ [ dy f + / 

Jy>0 Jy<0 J Uz>0 Jz<0. 



f dy± f dy\\f - f 

Jy>0 Jy<0 \ Uz>0 J z<0. 




(15) 



where the + , (— ) sign is relevant for pp , (pp) collisions. The forward-backward asymmetry 
is then given by 

do- jdM 



A FB (M) 



da+/dM 



(16) 



Figure [| displays our results for (a) the Drell-Yan invariant mass spectrum and (b) forward- 
backward asymmetry in the electron channel with and without scalar leptoquark exchange 
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at the Tevatron. Here, we have employed the present rapidity coverage of the CDF detector 
as used in their Drell-Yan analysis|^5[ |y maa; | < 1. We have assumed a scalar leptoquark mass 
of 200 GeV and Yukawa coupling strength of A/e = 1. In this figure, the SM is represented 
by the solid curve, and the cases with left-, right-handed leptoquark couplings to u-(d)quarks 
correspond to the dashed(dot-dashed), dotted(dot-dashed) curves. We see that the influence 
of leptoquark exchange on this process in minimal, even for these large values of the couplings! 
It is clear that at the present level of statistics, the Tevatron experiments are not sensitive 
to leptoquark exchange in Drell-Yan production. 

We next examine the level of sensitivity that will be achievable at the Main Injector 
with 2 fb _1 and \fs = 2TeV. Following Ref. f46fl , we enlarge the rapidity coverage to 



\Y\ < 2.5 and construct 21 invariant mass bins, corresponding to 

4 bins in steps of 10 GeV in the range 40 < M < 80 GeV , 
5 bins in steps of 4 GeV in the range 80 < M < 100 GeV , 

5 bins in steps of 20 GeV in the range 100 < M < 200 GeV , (17) 
5 bins in steps of 40 GeV in the range 200 < M < 400 GeV , 

2 bins in steps of 100 GeV in the range 400 < M < 600 GeV . 

The bin integrated cross section and asymmetry are then obtained for both the SM and 
for the case with 200 GeV scalar leptoquark exchange. We determine the statistical error 

on these quantities in each bin, which are taken to be 5N = y/~N and 5 A = J(l — A 2 )/N. 

The bin integrated results for the SM, along with the error associated with each bin, are 
displayed in Fig. ^. To evaluate what constraints may be placed on the leptoquark coupling 
we then perform a x 2 analysis according to the usual prescription, 

' Qi — Qf M \ 2 



bins 



5Q l 
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Figure 4: (a) The lepton pair invariant mass distribution and (b) forward-backward asym- 
metry in Drell-Yan production for the SM (solid curve) and with 200 GeV scalar leptoquark 
exchange, assuming A = 1, with left- (right-) handed couplings to u-quarks, corresponding to 
the dashed (dotted) curve, and left-(right-)handed couplings to d-quarks (dash-dotted curve 
for both cases). 
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where Qi represents each observable quantity. The resulting x 2 distribution, summing over 
both observables, is presented in Fig. |6] as a function of A/e for the various scenarios of 
left-, right-handed leptoquark couplings to u- and d-quarks. We see that the 95% C.L. 
bounds (corresponding to A% 2 = 3.842) on A, are quite weak and are inferior to the present 
restrictions from low-energy data. 

We briefly summarize this section by pointing out (i) present Tevatron data analyses 
easily allow for scalar leptoquarks in the 200-210 mass range, but exclude vector leptoquarks 
and leptogluons. (ii) If 200 GeV leptoquarks exist and are responsible for the event excess 
observed at HERA, the Main Injector will observe them in both the pair and single produc- 
tion mode, and can confirm the values of the mass and coupling observed at HERA, (iii) 
However, since these production mechanisms are QCD processes, hadron colliders can not 
provide any information on the electroweak properties of leptoquarks, or on the chiralities 
of their couplings. The same conclusions hold for the LHC. 

6 Constraints from Precision Electroweak Measurements 

With masses of only ~ 200 GeV, which is not far above the top quark mass, one may wonder 
if such light leptoquarks have any influence on Z-pole precision measurements. We first 
examine the oblique parameters |47|, [48], following the conventions of Peskin and Takeuchi, 
and in particular the extension by Maksymyk, Burgess and London to the enlarged parameter 
set S,T,U,V,W and X. The additional parameters V, W, and X need to be included when 
new light particles (with masses of order the electroweak scale or less) are introduced. Since 
we are assuming that the various leptoquark multiplets are degenerate, apart from higher 
order corrections, they do not contribute to the parameter T at one-loop (we recall that 
T is a measure of the mass splitting between the particles in a weak isospin multiplet). 
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Figure 5: Bin integrated lepton pair (a) invariant mass distribution and (b) forward-backward 
asymmetry for Drell-Yan production in the SM at the Main Injector. The vertical lines 
correspond to the expected statistical error in each bin. 
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Figure 6: y 2 fits to Drell-Yan production including the effects of a 200 GeV scalar leptoquark 
for each type of leptoquark coupling as labeled. The 95% C.L. constraints are obtained when 
X 2 = 3.842. 
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However, we do expect finite one-loop corrections to the other parameters, which can be 
straightforwardly calculated for each of the leptoquark multiplets. Here we present the 

results for the case of the F = type leptoquarks, R 2 l and R2l,r, in Fig- 0> which displays 
the shifts in the remaining oblique parameters as a function of the leptoquark mass. The 
results for the other leptoquark cases are found to be quite similar. As can be easily seen 
from the figure, leptoquarks in this mass range do not make appreciable contributions to the 
oblique parameters. 

In addition to oblique corrections, it is possible that relatively light leptoquarks can 
lead to substantial vertex corrections, e.g., in the case of Z — > e + e~, where first generation 
leptoquarks and quarks may contribute in a loop. This case has been previously examined 
by several authors^]. However, as shown by both Eboli et al. and Bhattacharyya et ai, 
the fact that first generation leptoquarks couple only to u- or ci-quarks leads to a substantial 
suppression of their potential contribution to this vertex. For leptoquark masses of order 
200 GeV, Yukawa couplings of order unity cannot be excluded by these considerations. We 
see that these constraints are much weaker than those imposed from low-energy data. 



7 e + e Colliders 

There are several ways in which leptoquarks may make their presence known in e + e~ 
collisions p)0"|, |5"T| , |52"| . At center of mass energies below the threshold for pair production, 



the existence of leptoquarks can lead to deviations |50|, ^2fl in the cross section and angular 
distributions for e + e~ — > qq. This may be particularly relevant when ^fs is comparable to 
the leptoquark mass as would be the case at LEP II if a 200 GeV leptoquark did exist. The 
origin of these modifications is due to the t— (u— )channel exchange of the F = 0(2) lepto- 
quark and is thus proportional in amplitude to the square of the unknown Yukawa coupling. 
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Figure 7: Shifts in the oblique parameters as functions of the leptoquark mass for the 
R 2 L and (b) R 2 l or R 2 r cases. 
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Including such terms results in the e + e~ — > qq tree-level differential cross section given in 
the Appendix. Note that for first generation leptoquarks either the uu or/and dd final state 
may be influenced, depending on the leptoquark species being exchanged. There will be no 
effect on ss, cc or bb final states. We now examine the sensitivity of the cross section to the 
value of the scaled coupling A. As an example, we consider the case of a 200 GeV leptoquark 
coupling to 4 with Xl = 1 at LEP II with \/~s = 190 GeV. Since d-quarks cannot easily 
be distinguished from any of the other light flavors, nor can quarks be differentiated from 
antiquarks without some difficulty [|53], we have symmetrized the expression in the Appendix 



with respect to cos 8 and summed over the possible light quark final states. Owing to this 
symmetrization, independent sensitivity to new t— versus u— channel exchange is lost and 
we can no longer distinguish F = from F = 2 type leptoquarks. Recall that in this sample 
case only di is assumed to couple to the leptoquark, hence this flavor summation will sig- 
nificantly degrade the sensitivity to the Yukawa couplings in this process. Figure |8| displays 
the resulting angular distributions for the SM and scalar (and for completeness, vector) lep- 
toquarks. From this it is clear that even large values of the Yukawa coupling do not lead 
to sizeable changes in the shape and size of the cross section. Of course there is nothing 
special about our sample case of leptoquarks coupling to di and we expect modifications of 
a similar size when the leptoquark couples instead to d^ or Ul,r- 

We now estimate the potential sensitivity of LEP II cross section measurements to 
non-zero values of the Yukawa coupling by generating a Monte Carlo data sample of qq 
events at 190 GeV with an integrated luminosity of 500pb _1 . The cross section is divided 
into 10 cos# bins of identical size and no additional cuts are applied. We then determine 
the sensitivity to the size of the Yukawa couplings assuming that only one of the quark final 
states couples to the leptoquark by performing a x 2 procedure. The results of this analysis 
are shown in Fig. |9] for leptoquarks of either spin and for both coupling helicities to up and 

31 




Figure 8: Symmetrized angular distribution of qq final states in e + e annihilation at 190 
GeV with z = cos 6. The solid curve is for the SM while the dotted(dashed) curve includes 

the contribution of a 200 GeV scalar (vector) leptoquark coupling to d L with A = 1. 
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down quarks. For scalar leptoquarks we expect that the LEP II limits on A will lie in the 
0.5-0.8 range which is somewhat less restrictive than those that already exist due to the APV 
data and the 7r — > ev decay. The expected bounds are also about a factor of 5 larger than 
the typical values necessary to explain the HERA excess in terms of leptoquarks. The cases 
where the leptoquarks are in a multiplet and can couple to both u and d quarks do not lead 
to any substantial improvement in these constraints when these contributions are combined. 
It thus appears that LEP II will be insensitive to any leptoquark consistent with the HERA 
data. 

At higher energy e + e~ colliders such as the NLC, leptoquark pairs can be produced 
directly and their properties examined ||50|| in detail. As we will see, this allows us to easily 



identify which leptoquark is being produced. In what follows we again limit our discussion 
to the spin-0 case since it is directly relevant for the HERA events. The pair production 
differential cross section is given by 

% = ^ 3 (1-^ 2 ){2E(« + «)W^ (19) 

S 2 



- f E + 4) + ~^r« - <)} C t P t + —(A! + A 4 R ) 



t , 



where the sum extends over 7 and Z exchange, (3 — yl — 4m 2 /s, C\ = Qlq, C 2 = 2[T 3 — 
x w Q]lq(V^G fM\ I 'Ana) 1 / 2 , and v , a , , and Pj are defined in the Appendix. 

When large Yukawa couplings are present, the exchange of u or d quarks in the 
t— channel can seriously modify the pair production angular distribution away from the con- 
ventional sin 2 9 dependence leading to an appreciable forward-backward asymmetry. How- 
ever, for A < 0.1 this asymmetry is found to be below the 0.2% level at a ^Js = 500 GeV 
collider when 200 GeV leptoquarks are being produced. Thus, in the limit that A ~ 0, the 
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Figure 9: x 2 fits to the SM angular distribution for e + e~ — > qq at 190 GeV including 
the effects of a 200 GeV leptoquark coupling to (a) u- or (b) d-quarks. In both cases the 
dotted(dashed) curve corresponds to a scalar leptoquark with a left(right)-handed coupling 
while the dash-dotted(solid) curve corresponds to the vector leptoquark case with left(right)- 
handed couplings. The 95% CL limits are obtained when x 2 — 3.842. 
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cross section for a particular final state, such as ££jj, depends solely on the electroweak 
quantum numbers of the members of the relevant leptoquark multiplet and their branch- 
ing fractions to charged leptons for a fixed leptoquark mass. Further information on the 
leptoquark electroweak quantum numbers may be obtained from examining the left-right 
polarization asymmetry, defined as usual as 

Table 5 summarizes the production cross sections for each final state as well as the polar- 
ization asymmetry associated with the ££jj final state for all scalar leptoquark multiplets. 
It is clear from the Table that by measuring the rates for each final state channel in ad- 
dition to the polarization asymmetry that the identity of the produced leptoquark would 
be straightforward to obtain assuming the design luminosity of 50 fb^ 1 . We note in passing 
that since the Yukawa couplings of the ~ 200 GeV LQ's are so small it is possible that LQ 
pair bound states(LQ-onium?) can form in the mass region near 400 GeV. These states can 
be produced in a number of ways, such as H^VT-fusion, and can only be explored in detail 
at a lepton collider such as the NLC although the lowest lying states can be produced at a 
hadron collider through gluon fusion. 



Single leptoquark production at the NLC is also possible via collisions |5^], with 
a production rate which is quite sensitive to the electric charge of the leptoquark. The 
amplitude for this process is proportional to the Yukawa coupling A and results in a cross 

section which is not significantly different in magnitude from that of pair production if A is 
not far from unity. If both electron and photon beam polarization is available, asymmetries 
can also be used to determine the leptoquark's quantum numbers as has been demonstrated 
by Doncheski and Godfrey [[54]]. Table 6 shows the production rate for a 200 GeV leptoquark 
in 7e collisions at the NLC for a luminosity of 50/6 -1 using either the backscattered laser 
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Leptoquark 


U33 




wjj 


ALR 


SlL 


1.88 


3.77 


1.88 


-0.618 


Sir 


7.53 


0.0 


0.0 


-0.618 


Sir 


120.4 


0.0 


0.0 


-0.618 


S3L 


192.2 


3.77 


1.88 


0.931 


R2L 


181.0 


0.0 


80.4 


0.196 


R2R 


261.4 


0.0 


0.0 


-0.141 


R2L 


47.6 


0.0 


33.2 


0.946 



Table 5: Cross sections for the three leptoquark pair decay channels in fb at a 500 GeV NLC 
assuming complete leptoquark multiplets with a common mass of 200 GeV. The polarization 

asymmetry in the ££jj channel is also given. In all cases A 1 is assumed. 

or Weisacker- Williams photon spectra. In obtaining these results, the hadronic content 
of the photon has been ignored; its inclusion would somewhat increase these rates. The 
backscattered laser approach has the advantage of a harder spectrum (although it cuts off at 
x ~ 0.84) and both beams can be polarized. In the Weisacker- Williams case an additional 

factor of 2 is included since both / ye ± collisions are possible. Even for this small value of A the 

production rates are at an observable level. It is clear that by using these rates together with 

the use of beam polarization the quantum numbers can be determined in a straightforward 
manner. 

8 Unification with Leptoquarks 

At this point one may wonder how scalar leptoquarks of the F = type would fit into a 
larger picture. As we saw earlier, both R2l,r can be embedded into a 45, 45 representation 

of SU(5) while R 2 l, which has less exotic electric charges, can be placed in a 10 or 15. 
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Leptoquark Backscattered Laser Weisacker- Williams 



S 1L 212. 56.8 0.5 

S 1R 212. 56.8 1 

S 1R 109. 24.4 1 

S 3L 430. 106. -0.75 

R 2L 332. 79.5 1 

R 2R 381. 92.6 1 

Ror 49.1 13.1 1 



Table 6: Rates for single leptoquark production in 7e collisions at a 500 GeV NLC assuming 
complete leptoquark multiplets with a common mass of 200 GeV and an integrated lumi- 
nosity of 50/6 -1 . In all cases A = 0.1 is assumed and a pr cut on the quark jet of 10 GeV 
has been applied. The charged lepton branching fraction for the produced multiplet is also 
given. 

In a SUSY extension, where one normally adds complete multiplets to automatically insure 
coupling constant unification, we would thus need to add either a 10+10 (15+15) or a 
45+45 at low energies. Here the barred representation is introduced to avoid anomalies and 
to guarantee that the fermionic components are vector-like with respect to the SM gauge 
group. As is well-known, the addition of extra matter representations delays unification and 
brings the GUT scale much closer to the string scale. A short analysis shows that adding 
complete 15+15's or 45+45's would lead to a dramatic loss of asymptotic freedom(AF) 
at one loop (i.e., Pi > 0) and, in the later case, both R2l,r would be present in the low 
energy spectrum. We would then need to explain why only one of the chiral couplings was 
present as well as the generational structure of the couplings by the imposition of some 
extra symmetries. The addition of the 10+10 to the usual MSSM particle content does 
not lead to either the loss of AF at one loop or to the problem of suppressing one of the 
chiral couplings. Interestingly, these general considerations tell us that the only leptoquarks 
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consistent with both SUSY and unification within standard SU(5) are R^l and Sil,r, the 
later being the familiar leptoquarks of E e string-inspired models|J. In both cases the QCD 
beta-function is found to vanish at one loop (in the E 6 case three 5+5's are present with a 
different leptoquark for each generation). Since we can safely add only a single 10+10 at low 
energies, a realistic model would still need to explain the hierarchy of generation dependent 
coupling strengths. We note that in the light 10 together with the leptoquark will be a 



Q — — 1 isosinglet bilepton|55| and a Q = 2/3, color triplet, isosinglet diquark. 



Although the leptoquarks in the 10+10 are found to be consistent with both unifica- 
tion and AF considerations, we still cannot identify them with the source of the HERA events 
due to the nature of their SU(5) coupling structure. To form a SU(5) singlet in the product 
5j5jl0fc, only the antisymmetric terms in the i,j can contribute. This would imply that 



the leptoquark must couple in an antisymmetric fashion with respect to the generations |56 
so that the phenomenologically required e + d-type coupling would be prohibited. Thus the 
requirements of AF, SUSY SU (5) unification, and the addition of complete multiplets do 
not simultaneously allow any F = leptoquarks with couplings to only a single generation. 

How do we circumvent this result? One possibility is to surrender the assumption 
of the addition of complete SU(5) representations at low energies. This certainly allows 
us more flexibility at the price of naturalness but still requires us to chose subsets of 
SU(3)c x SU(2)l x U(l)y representations from the 10, 15 or 45 which maintain AF and 
unification. Except for the rather bizarre choice of adding a (2,3)(l/6) from 15 and a 
(1, 1)(1) © (1,3)(— 2/3) from a 10 at low energy, a short analysis shows that no other solu- 
tions were found to exist. Here the notation refers to the (SU(3)c, SU(2) L )(Y/2) quantum 
numbers of the representation. Thus apart from this exotic choice we find that our con- 
straints are sufficiently strong as to disallow any F = leptoquarks at low energy in the 
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SUSY SU(5) context. 

It is clear that we must give up conventional £77(5) if we want a HERA-inspired 



leptoquark in a SUSY-GUT framework p5, |59|. Perhaps the most attractive scenario for 



this is the flipped £77(5) x U(l)x model|57]] wherein the SM fermion content is extended by 
the addition of the right-handed neutrino, z/ c , and the conventional roles of u c and e c are 
interchanged with those of d c and v c . Thus, u c lies in the 5, the d c lies in the 10 and e c is in 
an £77(5) singlet. In this case, completely different, and successful from the HERA point of 
view, leptoquark embeddings are now possible. For example, R2R can be placed in a 10+10 
without the difficulties associated with the cross generational couplings we encountered above 
since e c is an £77(5) singlet. In this case we would still need to impose some additional 
symmetries so that R2R could only couple to the first generation. Note that this 10+10 
would also contain the Sir leptoquark as well as an isosinglet bilepton with Q — 1. The 
Sir may also show up as a separate resonance in e~p collisions as discussed above if it has 
Yukawa couplings of order A ~ 0.1 once sufficient luminosity is accumulated. R2L can lie in 
a 10+10, but would require cross generational coupling as above since both L and u c are in 
the 5. On the otherhand, R2L now lies in a 45+45 and is excluded by the AF constraints. 
It thus appears that the flipped £77(5) x U(l)x scenario provides a natural embedding for 
at least one of these F = leptoquarks, R2R, and may predict the simultaneous existence 
of an F = 2 leptoquark. Other GUT groups may provide phenomenologically successful 
embeddings for the other F = leptoquarks. The extension of the spectrum to include the v c 
field may introduce some new additional interesting phenomenological implications for these 
leptoquarks since their interactions now extend beyond those described by the Lagrangian 
in Section 2 and may possibly yield excess events in the charged current channel. 

If LQ's really exist in a SUSY framework then their fermionic partners, leptoquarkinos(LQ- 
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inos) must also be observed [p0|. What are the masses of such states? If the LQ-ino is lighter 
than the LQ then the partner will decay directly to the LQ plus the LSP, since the LQ 
Yukawa is so small, unless the two states are nearly degenerate. This implies that the LQ- 
ino is the more massive of the two states, thereby decaying into the LQ plus LSP. Since the 
fermion pair cross section is larger than the scalar pair cross section at the Tevatron when 
both are color triplet states, it is clear that the mass of the LQ-ino must be substantially 
heavier than 200 GeV in order to reduce the anticipated number of eejj+ $t events to an 
acceptable level. 

What can we say about these leptoquarks in the non-SUSY context? Here we can 
only be more speculative. As is well-known, unification attempts without SUSY using only 
the SM particle content are doomed to failure in that they predict a too small value of the 
unification scale, implying a rapidly decaying proton, and lead to values of a s (Mz) which are 
smaller than the experimentally determined value by many standard deviations. One is led to 
consider the general question of whether one could add ad hoc sets of additional (non-SUSY) 
particles with masses at the electroweak scale to those which already exist within the SM to 
get unification at a higher scale and a proper value of a s (Mz) by sufficiently modifying the 
SM beta functions. A short consideration shows that this is difficult to arrange. At one-loop, 



the modified beta functions must satisfy the so-called "B-test" |61]: 



b — b 

B = ^ — ^ = 0.719 ±0.01 ±0.04, (21) 
b 2 - h 



where Bsusy = 5/7 = 0.714 clearly satisfies the test. In earlier work[^2j, many additional 
particles with a wide range of strong and electroweak quantum numbers were added in many 
thousands of combinations in order to attempt to satisfy these constraints with only two 
dozen candidates surviving (see Table 1 in Ref. ||62||). Given the higher precision of current 



data, at least several of these survivors could now be eliminated leaving a very short list. 
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A survey of this list shows that there is only one case with scalars which have the correct 
quantum numbers to be consistent with leptoquarks of any kind. Interestingly, this case 

corresponds to a pair of R 2 l leptoquarks with the Higgs sector of the SM augmented by 
an additional doublet. This scenario was first discovered in the analysis of Murayama and 
Yanagida|6S| and is quite unique with B = 0.693. Figure |Tl] shows a two-loop renormalization 



group(RGE) analysis of this particular case. Interestingly, for sin 6 W (MS) = 0.23165 and 
= 127.90 we obtain a s {Mz) = 0.123 and a proton lifetime [|64| of 10 32±1 yrs, which is 
close to the present limit in the e + 7r° modc|65||. 



As a final comment we note the oft-neglected problem of R-parity violation within 
a GUT context. A term in the superpotential of the form \ijk5i5jl0k would generate all 
of the usual lepton and baryon number violating R-parity violating terms simultaneously 
with comparable Yukawa couplings and would lead to a very rapid proton decay. It may 
be possible, however, to forbid such an R-parity violating term at the renormalizable level 
while having them arise as higher dimensional operators |)6). It is not clear, however, that 
an interaction of the type e c du can be generated in this approach with a Yukawa coupling 



of the right magnitude to explain the HERA events [57 



9 Conclusions 

In this paper we have considered the detailed phenomenological implications of interpret- 
ing the excess of events observed by both HI and ZEUS at HERA as the production of 
an s-channel scalar leptoquark resonance. First, we demonstrated that the DO leptoquark 
search data strongly indicate that this resonance could be neither a leptogluon nor a vec- 
tor leptoquark with a mass near 200 GeV due to their much larger pair production cross 
sections at the Tevatron. Secondly, we showed that the HERA data itself, in particular the 



41 




Figure 10: Two-loop RGE evolution of the model with the SM particle content together with 
a pair of R 2 l type leptoquarks and an additional Higgs doublet. 
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apparent lack of a signature in the e p channel even with the low accumulated luminosity, 
supports the idea that the leptoquark is of the F = type. We also showed that future 
HERA measurements in all e^ R p channels will allow a determination of the leptoquark's 
quantum numbers if sufficient luminosity and polarized beams become available. Thirdly, 
we analyzed the sensitivity of the Drell-Yan process at the Tevatron, dijet production at 
LEP 11 as well as precision electroweak measurements to the existence of leptoquarks. In all 

cases we found little sensitivity to leptoquarks with masses near 200 GeV with values of A 
near 0.1. Fourthly, we found that the single leptoquark production process at the Tevatron 
Main Injector may provide an independent determination of its Yukawa coupling, provided 
that sufficient integrated luminosity is obtainable, while pair production at the NLC allows 
one to directly determine all the leptoquark quantum numbers. Lastly, we saw that lepto- 
quarks can be embedded into a GUT structure both with and, surprisingly, without SUSY. 
Successful SUSY unification combined with the requirements of asymptotic freedom forced 
us to look beyond standard SU(5) to the flipped SU(5) x U(l)x model where R 2 r can be 
embedded. Without SUSY, a model with 2 Higgs doublets and a pair of R 2 l leptoquarks 
was found to unify near 10 15 GeV and led to reasonable values for both a s (Mz) and proton 
lifetime. We hope that this excess of events at HERA is confirmed by enlarged data samples. 
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10 Appendix: Relevant Formulae 

Here we collect the relevant formulae for the processes described in the text. 

• Vector Leptoquark Single Production at Hadron Colliders 

The k dependent parton-level cross section for the single production of vector lepto- 
quarks at hadron colliders is given by (z = cos 6) 




(22) 



with 




8 



- a{2x^XQ — x^m 2 ) + 2bx\X2 + 2c(xix 6 + X2X5 — 0:3X4) — 2dx±m 1 , 



v 2 = 



(2x 4 + m 2 ) 



2 



(23) 



^3 = 



£1(2x4 + m 2 ) ' 



t> 4 = 32x 4 x 5 + 16x 2 a;5 — I6X1X4 — 32xiX 6 + I6X3X4 — 16/txix 2 , 




and the definitions, 



a 



= -1 + 



2kxs 
m 2 



b=- 



2-2k + k 2 , 



(24) 



c = 



2 + Ak + 



(1 + K+ K 2 )X 3 
2 

wr 




(1 + k)x 3 
m 2 
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with 



Xi = 
%3 = 
%5 = 
t = 



is 



x 2 = 

= \(s - m 2 ) , 
Us - m 2 )(l + z) , u = -\{s - m 2 ){l - z) 



^m" -u) , 



(25) 



• SM and Leptoquark Contributions to Drell-Yan Production 



The SM contribution to the even and odd kinematic functions are 



K(sm) = 2(1 + z 2 ) + a\a\) {v\v{ + oJoj)P y , 



q{SM) 



1,3 



where the sum extends over the 7 and Z. Here we define 



s 2 [(s - M 2 )(s - M 2 ) + MjTjMjTj] 
ij ~ [($ - M 2 ) 2 + M 2 Y 2 ][{s - M 2 ) 2 + M 2 T 2 ] ' 



(26) 



(27) 



with M(r) being the masses(widths) of the gauge bosons, and the couplings are normalized 

as 



v} = Q f , 



a ' f = 0, 



V2G F M? 
Ana 



1 1/2 



(T 3 -2Q f x w ), a z f 



V2G F M 2 Z 
Ana 



1/2 



n. 



The scalar leptoquark contributions to these kinematic functions are 



F + 



- V 



^ \« + <M + a\)\ 2 Ltq + (vl - ai)(vi - a*)A^ 

i=f,Z 



i 2 



(28) 



+ 



u 



s(i-m 2 ) s(u-m 2 )_ 



45 



+- 



X L,q + X R,q 



i-2 



+ 



u 2 



(i-m 2 ) 2 {u-m 2 ) 



(29) 



7 q{LQ) = TE[(< + <)K + <)i + (<-«e)(<- fl ;)4 

i=j,Z 



s(u — m?) s{t — m 2 ) 



1 

±- 

4 



X L,q + X R,q 



u 2 



_{u-m 2 ) 2 (t-m 



2\2 



Here, 



s(s - Mi) 



(s - M 2 ) 2 + MfT 2 ' 



(30) 



and the top sign in the anti-symmetric function corresponds to the F — 0, S-type leptoquark 
exchange while the bottom sign is for the F = — 2, P-type leptoquark. Here, q represents 
either u- or d-quarks, depending on the coupling structure of the exchanged leptoquark. 

• SM and Leptoquark Contributions to e + e~ — > qq 



da 

dz 



37TQ; 2 



C 



Y, [i<4 + <ai)(vy q + 44) (i + z 2 ) + 2(viai + 4<)Ka? q + 4 a i>_ 



VI 



(31) 



x 



+ 



u 



s(t — m 2 ) s(u — m 2 ) 



± 



u 



s(u — m 2 ) s(t — m?) 



C 



X L,q + X R A 



U 



(t — m 2 ) 2 (u — m 2 ) 2 



± 



u 



(u — m 2 ) 2 (t — m 2 ) 2 



with the coupling normalizations and , Pj as defined in Eqs. (27,28,30), and the ± sign 
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is for 



± for < 



S, (F = -2) 
R, (F = 0). 



(32) 
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